{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 86,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "from scipy.integrate import odeint\n",
    "%matplotlib inline\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Polymer Cross-linkg reaction\n",
    " \n",
    " HPAM react with chromium ion as below:\n",
    " \\begin{equation}\n",
    " a[Cr]+b[Py]\\to c[Cr-Py]\n",
    " \\end{equation}\n",
    " \n",
    "Which is an irreversible reaction and the reaction rates are given by:\n",
    "\n",
    "\\begin{equation}\n",
    "   {{R}_{Cr}}=-{{k}_{f}}\\cdot {{{c}_{Cr}}^{1.22}}\\cdot {{{c}_{Py}}^{0.53}}\\cdot {{{c}_{H}}^{-0.9}}\n",
    "\\end{equation}\n",
    "\\begin{equation}\n",
    "  {{R}_{Py}}=-{{k}_{f}}\\cdot {{{c}_{Cr}}^{1.22}}\\cdot {{{c}_{Py}}^{0.53}}\\cdot {{{c}_{H}}^{-0.9}}\n",
    "\\end{equation}\n",
    "\\begin{equation}\n",
    "  {{R}_{Cr-Py}}={{k}_{f}}\\cdot {{{c}_{Cr}}^{1.22}}\\cdot {{{c}_{Py}}^{0.53}}\\cdot {{{c}_{H}}^{-0.9}}\n",
    "\\end{equation}\n",
    "\n",
    "Where Hydrogen ion (H+) concentration is calculated from measured PH value using (G. Paul, 2006):\n",
    " \\begin{equation}\n",
    "{{c}_{H}}={{10}^{-pH}} \n",
    "\\end{equation}"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Input Parameters\n",
    "Temperature 20 ℃:${{k}_{f}}=1.13\\times {{10}^{-3}}\\text{  }{{h}^{-1}}{{(mol/L)}^{0.15}}$<br />\n",
    "Temperature 35 ℃:${{k}_{f}}=3.16\\times {{10}^{-3}}\\text{  }{{h}^{-1}}{{(mol/L)}^{0.15}}$<br />\n",
    "Constant pH=4.9<br />\n",
    "Initial Concentration: ${{c}_{Cr}}=1\\text{ , }{{c}_{Py}}=60$"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 87,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "kf_20=1.13e-3\n",
    "kf_35=3.16e-3\n",
    "pH=4.9\n",
    "cCr_0=1\n",
    "cPy_0=3\n",
    "cCrPy_0=0"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 83,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "def Reaction_rate(C,t,k,pH):\n",
    "    cCr=C[0]\n",
    "    cPy=C[1]\n",
    "    cCrPy=C[2]\n",
    "    cH=10**(-pH)\n",
    "    \n",
    "    dCr_dt=-k*(cCr**1.22)*(cPy**0.53)*(cH**-0.9)\n",
    "    dPy_dt=-k*(cCr**1.22)*(cPy**0.53)*(cH**-0.9)\n",
    "    dCrPy_dt=k*(cCr**1.22)*(cPy**0.53)*(cH**-0.9)\n",
    "    return [dCr_dt,dPy_dt,dCrPy_dt]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 84,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAisAAAGkCAYAAADwoGW4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3Xl8VNX5x/HPE4IkhFUjIqIsioCC/lgU3FqQCkrBBdcg\nRcGKS8WK+9ZWweIOihataMWFTVErouIG7iugLQpVkcUVERFBkP35/XEmyWQySSaTCZPA9/163dfc\nOffcc58ZjHlyz7nnmLsjIiIiUlVlpDsAERERkdIoWREREZEqTcmKiIiIVGlKVkRERKRKU7IiIiIi\nVZqSFREREanSlKyIiIhIlaZkRURERKo0JSsiIiJSpSlZERERkSpNyYrIdsbMWprZP83sCzP71cx+\nNrM3zexCM8uqpGs2M7OtUdtmM1tqZk+a2YGVcc1kmFkjM7vNzBaY2Voz+8XMZpvZNWZWvwLt1oj5\n/FvM7Gsze97MDk/lZxDZEZnWBhLZfpjZ74HHgPXAw8DHwE7A4cCJwHh3P7cSrtsMWAxMBJ4DagBt\ngfMj1+/q7v9N9XXLw8wOisRWG3gUmBM51Bk4DXjb3Y9Osu0awCZgRqTtDKAl4fPnAke7+8sV+gAi\nOzAlKyLbCTNrDvwX+BI40t2XxxxvCfze3e+Kc64BO7n7hjjHzgT+5e4l3omNSlYudfdRUeV9gGnA\nP939vCQ+Vux1agO7ufvicp5XH/iEkET81t0/jzm+K3C2u48spY0sd19fwrH8ZOUOd784qvz/gLnA\ns+7etzwxi0ghdQOJbD+uAHKAs2ITFQB3X5SfqES6KsaYWX8z+5hwJ6ZXCe16ZEvGzMhri8h1F5vZ\nU7GVzKxWpLvqnjLa2xVYaGavmFmemdVKMI5zgd2BYbGJCoC7/xCdqJjZEjObZmY9zewDM/sVGJLg\ntaLb/QhYReHnf8vMZserG+m2e6a81xDZEShZEdl+9AEWuft7CdbvAYwCJgN/BpZUQkz7RF5XRF4f\nBY4xswYx9Y4F6gCPlNHed8AlhK6VR4HvzOyuyB2M0vQFfgWeSDBuB9oQurVeBC4EPkrw3AJmlgvU\nB36MFD0MdDCzfWPqHUJIaMr6/CI7JCUrItsBM6sL7AHMK8dp+xK6RG5w93tTNKaktpntEhnI+ltg\nPOEX/+OR4w8TxrCcEnPeAGCJu79dWuPuvtHd73D3A4GDgUlAf2Cumc0xs/NKGCjbFvjM3TeX47Ps\nDZzm7te4+zh3fz2Bc7Iinz/XzLpQ+Lkfi7xOATYSPm+0AcBq4OlyxCeyw1CyIrJ9qBd5XVOOc151\n909jC82sQeQX7i5mtgtQN1K+S8yWHafN64EfgGWELqAWwOXu/jRApAvmPeD0qOs1BI4m3ClJmLvP\ncfc/Ebp3+hPuXtxFuNvyiJk1japej/J9NwCLkxgUew7h8y8H3iEkVLe4+z8iMa8CpkfiBQrGu5wM\nPBFvzJCIQGa6AxCRlFgdea1bjnOWlFD+IdAsTvkPUftOSEyGx9S5j3A3YSthrMYn7r4pps7DwF1m\ntqe7f0W4y5JJVLJiZrvFnPNzSYNb3X0jMNnMHgfOA24jJAOPA19Hqq2mfN8NhAHDRUQSq52iita5\ne3QS9CRwD+H7WQ3Md/dfY5p5GOhnZl3d/V1CorYL6gISKZGSFZHtgLuvMbNvgXblOC32l2i+/kD0\nXZNewKXA7wCLKl8U59zP3X1mnPJok4HRhLsrN0VeZ8cMfP2O8AvfIq+DCL/kizGztpHjA4DdCE/9\n3A/Miqr2P+BAM8ssR1dQvO9nGnBYZN+BByg68ParBD7/84QxPAOAdyOv37j7qwnGJbLDUbIisv2Y\nDpxtZl3KMci2GHd/J/q9me0ZKZ8V/4xyt/+TmT0LnG5mEwm//C+Mqfa7mPefxMRUjzA3yiCgC6GL\nZzJwv7t/EOeyzwBdCXPNTKlA+H8GogcHf1PeBtx9s5lNBvLM7FrC4N9ij5OLSCGNWRHZftwCrAPu\nN7NGsQfNbG8zi00K0uURYH/gVmAzMQmEu8+M2b4HMLM6ZvYo4c7LPYTupsHA7u5+TgmJCsC9hHE0\nt5tZq9iDkQHB15QVtLvPjYmr2JifBD1CeKLpn4S7WBOSbEdkh6A7KyLbCXdfZGb9CXcYFphZ9Ay2\nhwEnAf9KY4jRniUMiD0ZeM7dV5RRP98uwFHAWMJEdQsSOcndV5nZCZHrfhRJePJnsO0I5AGlPomU\nSu4+28wWED7/f9394211bZHqSMmKyHbE3Z8xswOAywhzl5xLeFT2Y8K4k/vyq5L8RG8lXj7RNt19\nk5lNIQyIjTsWpQRfA3uU8xHk/Gu+b2btCN/N7wljRbYCnwE3U7QrJpnvp7znPAL8nfJ9fpEdkqbb\nF5G0MLNRhC6cxiU96bM9M7NLCAOM93L379Idj0hVpmRFRLa5yDT5XwHT3P2P6Y4nHczsE8LTQ0kt\nniiyI1E3kIhsM5EFA48ijJ/ZGRiT3oi2LTPLIXTP/Y4wnf/FpZ8hIqBkRUS2rf0Ik799DwxN0RT/\n1UljwpM/K4Hh7v5CmuMRqRbUDSQiIiJVmuZZERERkSpNyYqIiIhUaUpWRFLAzK4zs61mtnMK2trH\nzF40s1VmtsXMjk1FjCIi1ZWSFZHUSOUkaw8TpqK/GvgDMLsijZlZfTO7z8yWm9kvZjbTzDqU4/wm\nZvaYmf1kZj+b2b/NrEUJdY81szlm9quZLY0kcTVi6uSY2fVm9ryZ/RhJ8gaW0N5BZjbWzGab2UYz\n21JCvSwze8DM5kWSvDVm9pGZXWhmlf4gQVnJqpktMbNplXDdwyPX3ZJoohxZdmGqma00s7Vm9oaZ\ndUt1bCKppKeBRKoQM8siLLg3wt3HpqA9A54D2hPWDvoROB941cw6uvsXZZyfA7wK1AVuIKzjc3Hk\n/P9z95+i6h4DPAXMBC6IXPNaYFfgT1HN5gJ/AZYCHwHdSgmhN2HiuP8CXwD7llAvG2hLmE5/CWFm\n2kMJqzsfTJittjKVlaym/EmGyL/tXcAvQE6C5zQlrPS8iTBr7zrCYpAvmtmR7v5mquMUSQUlKyJV\nS/4ChD+nqL2TgUOAE939KQAze5wwxfz1lP1L/E/A3sBB7j43cv4MwvT9lxCSkXy3EZKPXu6+NVJ3\nDXCVmd3p7p9F6n1LmLV2uZl1AkpafBDCGkA3ufsGM7uLEpKVSNJ0aEzxfWa2GviTmV3s7svL+KzV\nzTnAHsD9hNWgE3EVUA/Y390XApjZ/cD/CIndQZUQp0iFqRtIpJKYWTMzW2hm/41MhlZW/b8R7go4\ncFvk9v6iCoZxIrAsP1EBiCwa+BhwnJnVTOD8D/ITlcj5nwKvAKdExd6WcGfjvvxEJWIs4f8zJ0Wd\nvynRxMHdf3D3DYnULcHSyGuD8p4Y+ffbamYXm9lFka6cdWb2qpntX4GYKszMGgIjCHeoypPYHg58\nmJ+oALj7r8A0oKOZ7Z3SQEVSRHdWRCpB5H/6M4EfgKOiu0tK8QTwE3AHMJHQffNLpL1MoH6Cl1/p\nhRModQDmxqnzPnA24U7FJyV8BgMOAB4o4fyjzCzH3ddGruMUrmQMgLt/Z2ZfR45XukjyVY/QLXQQ\n4e7PEmBhKaeV5QygDnA3kEW4i/GKmbV39x9i6u4S+d6KhEWcPwzNrAFQI7Y8jnWRhCLaDcB3hIUp\n/5pAG/lqESakK3aNyGsnQnebSJWiZEUkxcysDfAyYe2bo909ob983f3jSLfJHcBcd58YdfgwYFYi\nzQAtgC8j73cHXotTL3/hvCaUkKwQpsOvFVW3pPM/j1yHUuo2KTXq1OkHTIp6/wEwOOZuT3ntDezj\n7ssAzOwF4D3gCsJK1vkM+LSENhz4T0zZh0CzMq7thO664QUXCatqDyH8t+XFc6NSfQocHpVk5jsi\n8rpHeRoT2VaUrIikVntgCmFMSG93/yVF7X5EWE8mEcui9rOBeN0o6wm/XLNLaSf/WEnnR9cpq27d\nUq6TSjMJ31MDoAdwIOGuSEU8lZ+oALj7B2b2HmHwb3Sy4oRkaU2cNibEKetP6d9/vtiuwDHAs+7+\nSgLnxroH6As8ZmbXAGsJ45I6RY4nEo/INqdkRSR1DHiGkCwc7e7ryqifsMjdmZlJnPor4e5IrCzC\nL9fY7oXYcynl/Og6ZdUt7TopE+mWyf+enjSzq4CXzGyfCgywjdeF9Blh8HKsN9y9WDeLma2PLXP3\nd8obiJmdSnhaLKkxM+4+w8wuAG4idNkZ4c7Y1cCtRLodRaoaDbAVSR0HphK6DVL6qKyZ1TSz3RLc\non+uv6OwiyZaftm3pVx2JeFOSSLnfxdTHlu3tOtUpqmEOyvHpen6JTKz3AT/PaMfS74FeBzYHBkA\n3AxoGDm2l5nF+/6LiDwSvxvh6alOhNWfVxP++/2slFNF0kZ3VkRS6zJgCzDWzFa7++QUtXsoyY1Z\n+YjwBEisroRBlSX+coqMh5gHdI5zuAuwKGrcw0eEv9I7EzWJXeSXZ1Pg3gRirwz53RqJDk6Op1Wc\nsn0JA3cr4gPKP2ZlT0L30elx6s4l/Dt0LOvCkQG77+W/N7OjCHe/3iozapE0ULIiklpOGPxYF3jY\nzH5x9+kpaDfZMStTgRPNrJ+7PwnhL3rCo8TT3H1TfkUz2xOoHXk0Ofr8GyMTyOXPs9IaOJLwVz4A\n7j7fzP4HDDGzf0Y9jXQ+YYK2J8rxWcvNzHZx9x/jHDqb8G9SkVmAjzezJu7+beRaBxOStVEVaBOS\nG7NyfJzjeYTHyP8AfJNfaGb1CHe1vnP31SU1bmaHAicA/3D3eONtRNJOyYpIikXuSAwA/g08bma9\n3T2RuyKltZnsmJWpwEXAg5G5QVYQEogM4LqYuo8Av6Fo9/BYwi/858zsNsIMtsMI3T6xv6wvA54m\njBGZTBhs/CdgXEwChJn9iTAINv/pk2MjyRLAmPxfmma2F+GXMETu8EQGhgIsdfdHI/sDzOxcwne+\niJAs9iIkeNPc/dWSv6IyLQTeNLN7KHx0+QfCGI+kJTNmxd2LTdlvhUsnzIgZL3MC8CBwJmEJh/zv\n8zHCvCrLgHaEyeU+Aq5BpKpy9yq1AfsRfpi+IIxU/4Hw6GWfBM+vT5h7YDlhsNhMoEO6P5e27XsD\n/kbo/tk5qiwr8t/fz4QZYBNpp1mknWEpjC36Z2INYUK3Yj8ThG6mzXHKmxCecPop8ln+DbQs4VrH\nEgZuriNMyHYdUCNOvcWRzxlv2yuq3m8Jd2bi1ZsZVa8TMDnS7jrCGIwPgAuBjCS/t2aRa19MSPiW\nRNqeBbQr698/5vgi4Olt9d9epPyMSPnAqLIGwJOEOzC/EhKxvwM56fz50aatrM3cU75kRYVE1hcZ\nCrxDGJRXmzCL5m+AIe5+fynnGvAmxddB2Qsocx0UEZF8kcGri4FL3b2iXT4iUgFVLlmJJ5KEzAVq\nuft+pdQ7hfDXVfQ6KLmEQYTPuXtlL2YmItsJJSsiVUe1GLPi7m5mXxH/qYRocddBMbPHgNPNrKZH\nDSgU2VYsrKZc1hMpK/Xfp4hIcVV2nhUzq21mu5hZSzMbBhxDmMK8NKWtg1KbkpeXF6lspxIGpZa0\nfUtYHVmqFo9sIpJGVfnOyu2EUepQ+Ojj0DLOqcg6KCKVaQZlP3ocu3aMpJG7LyWxhQZFpJJV5WRl\nNGGmxiaEOQRqEH8q72gVWQdFpNK4+/fA9+mOQ0SkOqqyyYq7f0bh7JqPRlY6nUaYebMkSa+DYma7\nEOZlWELhIm0iIiJStiygOfCCx5+gsUKqbLISx1TgXjNr5e6fl1CnIuug9CL+yqgiIiKSmNOBialu\ntDolK4ms8ZH0OihE1vnYY49HmTatbTLxSRKGDRvG6NGj0x3GDkXf+ban73zb03e+bS1YsIABAwZA\nxdfMiqvKJStmtquHZd6jyzIJszH+CsyPlDUmJC4L3X1LpGrC66DEsR7g22/bsu++HalTJ5WfSkpS\nv359OnYsc901SSF959uevvNtT9952lTKMIoql6wA/4wswPU6YUroxoTbSq2Bi919XaTeTcBAQh9Z\n/gqz5VkHJS53+PBDOOKIlHwWERERqaCqOM/KZMJ6FucSFlEbBnwFHOvud0bVc8IjzYUF7lsJ87FM\nITzmfAthPZTupYxzKWZ2RdZnFRERkZSqcndW3P0xwkKGZdUbBAyKU/4zMCSyJeWDD5I9U0RERFKt\nKt5ZSTslK9tOXl5eukPY4eg73/b0nW97+s63L9ViIcNtwcw6AnPC6vYdWbkSGjZMd1QiUhm+/PJL\nVqxYke4wRKqV3Nxc9tprr7jH5s6dS6dOnQA6uXu8ZW8qpMp1A1UVc+bA78qaHF1Eqp0vv/yStm3b\nsm7durIri0iB2rVrs2DBghITlsqkZKUEs2crWRHZHq1YsYJ169bx6KOP0rat5lQSSUT+PCorVqxQ\nslKVaNyKyPatbdu2modDpJrQANsY2ZF5cpWsiIiIVA1KVmLk3xX+6iv4XmvkioiIpJ2SlRj771+4\nr8nhRERE0k/JSoz99ivcV7IiIiKSfkpWYkQnKxq3IiIikn5KVmLssQfsvHPY/+CDsLChiEh1s2jR\nIs455xz23ntvsrOzqV+/Pocffjhjxoxh/fpKWRhXpNLo0eUYZtC5M7z4IixfDl9/DXvume6oREQS\n9+yzz3LKKaeQlZXFwIEDadeuHRs3buTNN9/k8ssvZ/78+dx7773pDlMkYUpW4jjooJCsQLi7omRF\nRKqLJUuWkJeXR4sWLZg5cyaNGjUqOHbeeecxYsQInn322bjnujsbN26kVq1a2ypckYSoGyiOzp0L\n9zXIVkSqk5tvvpm1a9fywAMPFElU8rVs2ZKhQ4cCkJGRwYUXXsjEiRNp164dWVlZvPDCC9s6ZJEy\n6c5KHAcdVLivQbYiUp1Mnz6dli1b0qVLl4Tqv/LKKzz22GNccMEF5Obm0rx588oNUCQJSlbiaNIE\nGjeGZcvCnRX3MJZFRKQqW7NmDd988w3HH398wud89tlnfPzxx7Ru3boSIxOpGCUrcZjBwQfDtGmw\nahXMmwcHHJDuqEQkLTp3Dn+5VKbGjVPS57x69WoA6tatm/A53bp1U6IiVZ6SlRL87nchWQF4/nkl\nKyI7rGXL4Jtv0h1FQurVqweEOyyJUrePVAdKVkpwzDGF+88/D1dckb5YRCSNGjeuNteoW7cuTZo0\n4eOPP074nOz81VtFqjAlKyXYZ5+wLVwIb70Fq1dD5I8WEdmRVLNHAvv06cO4ceN47733Eh5kK1LV\n6dHlUuTfXdm8GV5+Ob2xiIgk4vLLL6d27dr88Y9/ZPny5cWOf/HFF4wZMyYNkYkkT8lKKWK7gkRE\nqrqWLVsyceJEFi1aRNu2bRk2bBgPPPAA99xzDwMGDGD//ffnf//7X7rDFCkXdQOVols3yMqC9etD\nsqJHmEWkOujbty///e9/ufXWW5k2bRr33nsvO+20E+3ateO2225jyJAhAJgZpv+pSTWgZKUU2dkh\nYZkxIzwM8PHH0L59uqMSESnb3nvvXeb6P1u2bNlG0YhUjLqByqCuIBERkfRSslIGJSsiIiLppWSl\nDK1awd57h/033wyPMIuIiMi2o2QlAdGPML/ySnpjERER2dEoWUmAuoJERETSR8lKArp1g1q1wn7+\nI8wiIiKybShZSUDt2iFhAfj6a/jkk7SGIyIiskNRspIgdQWJiIikh5KVBEUnK9OmpS8OERGRHY2S\nlQS1agVt24b9N98MqzGLiIhI5VOykiAzGDy48P2DD6YvFhERkR2JkpVy+MMfoEaNsD9+PGhZDRHZ\nni1dupSMjAwefvjhdIciOzglK+Ww227Qp0/Y//ZbeOGF9MYjIhLroYceIiMjo2DLzs6mdevWDB06\nlOXLl6c7PJGkVLlkxcw6m9ndZvaxmf1iZkvNbIqZtUrg3DPMbGucbYuZNUpFfNFdQf/6VypaFBFJ\nLTPjhhtu4NFHH+Uf//gHhx12GPfccw+HHnoo69evT3d4IuWWme4A4rgCOBR4HPgv0BgYCsw1sy7u\nPr+M8x34C7AkpnxVKoLr3RsaN4Zly8JTQT/8ALvumoqWRURS5+ijj6Zjx44ADB48mJ133pnRo0fz\n9NNPc+qpp6Y5usq1YcMGdtppJ8ws3aFIilS5OyvA7UAzd7/I3f/l7iOBIwiJ1ZUJtjHD3SfGbBtT\nEVxmJgwcGPY3bYJHH01FqyIilevII4/E3Vm8eDEAixcv5uSTT2aXXXYhJyeHQw45hOeee67UNsaP\nH09GRgb/+c9/ih0bOXIkmZmZfPfddwB069aNAw44gHnz5tGtWzdycnJo1aoVTzzxBACvvfYaXbt2\npXbt2rRp04ZX4iy89u233zJ48GAaN25MVlYW7dq148GYpxtee+01MjIymDJlCtdeey1NmzYlJyeH\nNWvWJPU9SdVU5ZIVd3/X3TfHlC0EPgHaJtqOmdUxs0r5fIMGFe4/8ICm3xeRqm9hZL6F3Nxcli9f\nziGHHMJLL73EBRdcwMiRI9mwYQPHHnssTz/9dIltnHTSSWRnZzNhwoRixyZOnMiRRx7J7rvvDoSu\nqJUrV9K3b1+6du3KrbfeSlZWFnl5eTz22GPk5eXRp08fbr75ZtauXcvJJ5/M2rVrC9pbvnw5Xbp0\nYebMmVx44YWMGTOGVq1acdZZZzFmzJhi1x8xYgTPP/88l112GSNHjmSnnXaq6FcmVYm7V4sN+Ap4\nvow6ZwBbgdWR1/XA08A+CbTfEfA5c+Z4Ig491D2kKe7vvZfQKSJSBcyZM8fL87Ne3YwfP94zMjJ8\n5syZvmLFCv/666998uTJnpub63Xq1PFvv/3WL7roIs/IyPC333674LxffvnFW7Zs6S1btiwoW7Jk\niZuZP/TQQwVl/fv396ZNmxa55ty5c93M/OGHHy4o69atm2dkZPiUKVMKyj799FM3M8/MzPQPPvig\noPzFF18sdp2zzjrL99hjD//pp5+KXCsvL88bNmzo69evd3f3V1991c3M99lnH9+wYUOyX5uUoayf\nm/zjQEevhBygKo5ZKcbMBgB7ANeWUXUd8CAwi5CwdAIuAd4ys47u/k2qYho8GN5+O+z/619w8MGp\nallEqpLOncMYtcrUuDHMnp269tydHj16FLw3M5o3b86kSZPYfffdef755zn44IM55JBDCurk5OQw\nZMgQrr76aubPn89+++0Xt+2BAwcyefJkZs2aRffu3QGYMGECtWvXpl+/fkXq1qlTh1NOOaXg/b77\n7kuDBg1o2rQpnTt3Lijv0qULAIsWLSooe/LJJzn11FPZsmULP/74Y0F5z549mTJlCnPnzi0S/5ln\nnqm7KduxKp+smFkb4G7gLaDUh/3d/XHCwNx808zsReB14Brg/FTFdcop8Oc/w9q1MGkSjBoVFjwU\nke3LsmXwTcr+zNk2zIyxY8fSqlUrMjMz2W233WjdunXB8aVLl9K1a9di57WNTNO9dOnSEpOVo446\nisaNGzNhwgS6d++OuzN58mSOP/54cnJyitRt2rRpsfPr16/PnnvuWaSsXr16APz0008A/PDDD6xa\ntYr77ruPf/7zn3E/X+xj2M2bN48br2wfqnSyYma7Ac8CPwEnu5d/dIi7v2Vm7wG/S2VsdeuGhOXB\nB2H1apgypehYFhHZPjRuXD2vcdBBBxU8DZRKGRkZ9O/fn/vvv5+xY8fyxhtv8O233zJgwIBidWvk\nz6KZYHn+/+K3bt0KwIABAzjjjDPi1j3ggAOKvM/Ozk74M0j1U2WTFTOrB8wA6gGHu3tFbsR+Beyb\nSMVhw4ZRv379ImV5eXnk5eUVq3v22YXT7o8cGWa4zayy36iIJCOV3TNVRbNmzfj000+LlS9YsKDg\neGkGDhzIqFGjeOaZZ3juuedo1KgRPXv2TFl8u+66K3Xr1mXLli0ceeSRKWtXUmPSpElMmjSpSNnP\nP/9cqdeskr9azawWMB3YB+jh7sV/qsqnJfBDIhVHjx6d8F8jhxwC3bvDrFlhYcMJE6CEPwJERKqM\n3r17c+edd/Lee+8VjBdZu3Yt9913Hy1atCixCyhf+/btad++PePGjePdd99l0KBBZGSk7uHLjIwM\nTjzxRCZNmsRVV13F/vvvX+T4ihUryM3NTdn1pHzi/QE/d+5cOnXqVGnXrHLJSuRx48eALsCx7v5+\nCfUaA/WBhe6+JVKW6+4rYur1Jgy0vaMy4r3++pCsAIwYAaefrrsrIpJeZfWYX3nllUyaNImjjz6a\nCy+8kJ133pnx48ezdOlSnnzyyYSuMXDgQC699FLMjNNPPz0VYRdx00038eqrr9KlSxfOPvts9ttv\nP1auXMmcOXOYOXMmK1asKLsR2W5UuXlWgFFAX+B5INfMTo/eourdBCwgPCWU7+3I1PyXmdkQM/sn\n8G9gKXBjZQR7xBGQP+j+iy/gkUcq4yoiIokra+bWRo0a8c4779CzZ0/uvvturr76arKyspg+fTrH\nHntsQm2dfvrp1KhRg9atWxd5sqesc80sofJGjRrx/vvvM3jwYJ566imGDh3KmDFjWLVqFbfccku5\nPq9Uf1XxHsCBhGe1+0a2WPmzETlhLpVok4HfA0cBtYHvgH8Cw909oW6gZFx/PeRPvjhiBAwYADVr\nVtbVRERKdsYZZ5Q4KDVa8+bNmTJlSql1mjVrxpYSlpfPzMzEzOIOrAWYlX/LOUb048nR4l0nNzeX\nMWPGxJ0ELt9vf/vbEmOU7UeVu7Pi7t3dvUZJW1S9Qe6e6e5fRpX91d07ufvO7p7l7i3cfWhlJioA\nhx0GRx0V9hcvBq2mLiLbuwcffJCtW7eWmKyIpFKVS1aqq+uvL9y/4QbYmJKViEREqpZZs2Zx9913\nM3LkSE444QT22muvdIckOwAlKylyyCHQq1fYX7IEHnooreGIiFSK4cOHc+mll9KxY8dSu2dEUknJ\nSgrF3l1Zvz59sYiIVIZZs2axfv16Xn755YJFC0Uqm5KVFOrSBY45Jux/+WUYbCsiIiIVo2QlxW69\ntfBJoJtvho8+Sm88IiIi1Z2SlRTbf3+45pqwv2VLWJ158+b0xiQiIlKdKVmpBFddBe3ahf0PP4Tb\nb09vPCK8Yfx5AAAgAElEQVQiItWZkpVKsNNO8MADkL9Uxt/+Bp99lt6YREREqislK5Xk4INh2LCw\nv2EDnHUWbI2db1dERETKVO5kxcxqm9lxZnaLmU03s3fN7B0zeyZSdpyZ5VRGsNXN8OGw995h/803\nYezY9MYjIiJSHSWcrJhZezMbDywDngL+BOxDWKPHgH2BCyLHlpnZeDNrn/KIq5HatWHcuML3l1wC\n776bvnhERKqT6667joyMor+mmjdvzuDBg5Nqr3nz5sUWaqwMFYlR4ksoWTGzKcCHQBvgOsJig/Xc\nvY27H+LuXd29NVA3cuw6oDXwoZlNqozAq4vu3WHo0LC/cSOccAJ8/XV6YxKR7duiRYs455xz2Hvv\nvcnOzqZ+/focfvjhjBkzhvXVaLbKeCs0Z2RkJL3K8rZanbk819mwYQOjR4+ma9euNGjQgOzsbFq3\nbs3QoUP5/PPPE2pj0KBBZGRkFGz169fn//7v/xg1ahQbt5O1XxJddXkr0NndS501xN23APMi2+1m\n9n/AFRULsfq7/XaYNw9efRWWLYPjj4fXXw93XkREUunZZ5/llFNOISsri4EDB9KuXTs2btzIm2++\nyeWXX878+fO599570x1m0j799NNid1uqqx9//JFevXrx4Ycf0qdPH04//XTq1KnDp59+yuTJkxk3\nblzCyWVWVhYPPPAA7s6qVat44oknuPTSS5k9ezYTJ06s5E9S+RJKVtw9L5nGI8lNUuduT2rWhMcf\nD4NuFy+GOXPCgNuJE2EbJfoisgNYsmQJeXl5tGjRgpkzZ9KoUaOCY+eddx4jRozg2WefLfF8d2fj\nxo3UqlVrW4SblJr5s25uB8444wz+85//8MQTT3D88ccXOTZixAiuyZ+0qwTr1q2jduSv3szMTPLy\nCn/dnnfeeXTp0oUpU6YwatQoGjdunPoPsA1VSnpqZr3N7L7KaLu6ys2FadOgTp3wfvJkuPHG9MYk\nItuXm2++mbVr1/LAAw8USVTytWzZkqH5/dKELpULL7yQiRMn0q5dO7KysnjhhReA8IvwkksuYa+9\n9iIrK4s2bdpwe5xJo1566SWOOOIIGjZsSN26dWnTpk2xX7J33XUX7dq1Iycnh5133pmDDjqIyZMn\nJ/UZY8eDPPTQQ2RkZPD2229z8cUX06hRI+rUqUO/fv348ccfy2zvoYceombNmlxxRWEngLtzxx13\n0K5dO7Kzs2ncuDHnnnsuq1atKnb+DTfcwJ577klOTg49evRg/vz5CX2O999/n+eee44//vGPxRIV\nCEnZLbfcUvD+zDPPpG7duixatIjevXtTr149BgwYUGL7Zka3bt2AkMSOHz+ejIwM/vOf/xSrO3Lk\nSDIzM/nuu+8Sij0dKuteWgfgrEpqu9pq1w4efbTwbso118CkHXpEj4ik0vTp02nZsiVdunRJ+JxX\nXnmFiy++mNNOO40777yT5s2bA9C3b1/uvPNOevfuzejRo2nTpg2XXXYZl1xyScG58+fPp2/fvmza\ntIkRI0YwatQojjvuON5+++2COuPGjePPf/4z7dq1484772T48OF06NCB9957L6nPWNJ4kKFDhzJv\n3jyuu+46zj//fJ555hkuuOCCUtu67777GDx4MFdffTU333xzQfmQIUO44oorOOKIIxgzZgyDBw9m\nwoQJHH300WzZsqWg3l/+8hf++te/0qFDB2677TZatmxJz549WbduXZmfY9q0aZhZqQlH7OfevHkz\nvXr1onHjxtx+++2ceOKJpZ6zcOFCAHbZZRdOOukksrOzmTBhQrF6EydO5Mgjj6zaC1O6e8o34Bpg\nS2W0XVkb0BHwOXPmeGX7+9/dIWwZGe4PP1zplxSRiDlz5vi2+lnfllavXu1m5ieccELC55iZZ2Zm\n+v/+978i5f/+97/dzPzGG28sUn7yySd7jRo1fNGiRe7ufscdd3hGRoavXLmyxGscf/zx3r59+3J8\nkkLXXXedZ2RkFClr3ry5Dxo0qOD9+PHj3cy8V69eRepdfPHFXrNmTV+9enWRc/v27evu7nfeeadn\nZGT4yJEji5z3xhtvuJn55MmTi5S/+OKLbmY+adIkd3f/4YcfvFatWn7ssccWqXfNNde4mRWJMZ5+\n/fp5RkaG//zzz6XWy3fmmWd6RkaGX3PNNXGP1a1b11esWOErVqzwL774wkeOHOkZGRneoUOHgnr9\n+/f3pk2bFjl37ty5bmb+cBm/iMr6uck/DnT0SvgdnegAW0mhq66CJUvCY81bt8IZZ8CmTWEdIRGp\nWjrf15llvyyr1Gs0rtOY2UNmV6iN1atXA1C3bt1yndetWzdat25dpOz5558nMzOzSJcRwCWXXMLU\nqVN5/vnnOf/882nQoAEATz31FIMGDYp716NBgwZ8/fXXzJ49m86dO5crtkSZGUOGDClSdsQRR3DH\nHXewdOlS2uWvfxJx6623csUVV3Dbbbdx8cUXFzk2depUGjRoQI8ePYp0I3Xo0IE6deowa9YsTjvt\nNF566SU2bdpU7Du66KKLGDlyZJkxJ/vvde6558Yt/+WXX9h1110L3psZhx12GA8//HBB2cCBA5k8\neTKzZs2ie/fuAEyYMIHatWvTr1+/csWxrSlZSQMzuPfeMPB27Nhwj+Wss8JMt+edl+7oRCTasl+W\n8c2ab9IdRpnq1asHwJo1a8p1Xn63T7SlS5fSpEkTcnKKzu/Ztm3bguMAp556Kg888ABnn302V155\nJT169KBfv36cdNJJBYnLFVdcwSuvvMLBBx/MPvvsQ8+ePenfvz+HHnooAJs2bWLlypVFrtOoUaNy\nP2a85557FnnfsGFDAH766aci5a+++irTp0/nyiuvLJaoAHz++eesWrUq7pgfM2P58uUAfPnllwDs\ns88+Rerk5uYWXLs00f9e+ftlyczMpGnTpnGPZWdnM336dNydWrVq0aJFC5o0aVKkzlFHHUXjxo2Z\nMGEC3bt3x92ZPHkyxx9/fLF/66pGyUqaZGTA3XdDrVowenQoO//8kLBcdFF6YxORQo3rVP5TFKm4\nRt26dWnSpAkff/xxuc7Lzs5O+ppZWVm8/vrrzJo1i2effZYZM2YwZcoUevTowYsvvoiZ0aZNGz79\n9FOmT5/OjBkzePLJJxk7dix/+9vf+Nvf/sbbb79N9+7dMTPcHTNj8eLF7LXXXuWKpUaNGnHLPXTz\nF2jXrh2rVq3ikUceYciQIcWSta1bt7LbbrsxceLEYucCRe5eVESbNm0AmDdvHocddlhC55T2lFaN\nGjUK7paUJCMjg/79+3P//fczduxY3njjDb799tuEx82kU8LJiplNK0e7+5RdRczCHCy1asFNN4Wy\nYcPCood33BEWRBSR9Kpo98y21KdPH8aNG8d7771XrkG2sZo1a8Yrr7zC2rVri/zFvWDBgoLj0bp3\n70737t257bbbuPHGG7n22muZNWsWRx55JBASopNPPpmTTz6ZzZs3c8IJJ/D3v/+dq666igMPPJCX\nX365SHuV+Zhtbm4uU6dO5bDDDqNHjx689dZbRa63995788orr3DooYeWmhzkfweff/55kYRnxYoV\nxe7mxNO3b19uvPFGHn300YSTlVQYOHAgo0aN4plnnuG5556jUaNG9OzZc5tdP1nleRroAKB9gls2\n8GVKI91OmcHIkWFl5nz33AO//S18U/XvPItIFXL55ZdTu3Zt/vjHPxZ0V0T74osvGDNmTJnt9O7d\nm82bN3P33XcXKR89ejQZGRkcc8wxQPEuFoADDzwQd2fDhg0Axbp4MjMzadu2Le7Opk2baNCgAUce\neWSRbadK/kutSZMmvPzyy/z6668cddRRRT7HKaecwubNmxk+fHix87Zs2cLPP/8MwO9+9zsyMzO5\n6667itQZnX+rvAxdu3bl6KOP5v777+fpp58udnzjxo1cdtll5flYCWnfvj3t27dn3LhxPPHEE+Tl\n5VWLSfYSvrPi7s0rMY4dmhlcdx00bw7nnhu6gt59Fzp1gsceg9/8Jt0Rikh10LJlSyZOnMhpp51G\n27Zti8xg+9ZbbzF16lQGDRpUZjt9+/ale/fuXHPNNSxevJgDDzyQF154gWeeeYZhw4bRokULAIYP\nH87rr7/O73//e5o1a8b333/PPffcw1577cXhhx8OQM+ePWncuDGHHXYYu+22G/Pnz+cf//gHffr0\nSdk4iXjdNaWVQ7iD8uKLL9KtWzd69uzJzJkzqVu3Lr/5zW8455xzuOmmm/joo4/o2bMnNWvW5LPP\nPmPq1KmMGTOGfv36kZuby6WXXspNN91Enz596N27Nx9++CEzZsxIuKvo4YcfplevXpx44on06dOH\nHj16kJOTw+eff87kyZNZtmwZt956a1LfSWkGDhzIpZdeiplx+umnp7z9SpHoY0PASUDdyngkqSps\nbMNHl0sze7Z7s2aFjzbXqOF+883umzalNSyR7cb2+uhytIULF/o555zjLVu29KysLK9Xr54feuih\nftddd/mGDRsK6mVkZPiFF14Yt421a9f6JZdc4k2bNvVatWp569atfdSoUUXqzJo1y0844QRv2rSp\nZ2VledOmTX3AgAG+cOHCgjrjxo3zbt26+a677urZ2dneqlUrv/LKK33NmjVlfo7rrrvOa9SoUaSs\nRYsWPnjw4IL348eP94yMjGL/nq+++qpnZGT4a6+9VuTc2EeN33//fa9fv75369bN169fX1B+//33\n+0EHHeQ5OTlev359P/DAA/2qq67yZcuWFTl/xIgRvscee3hOTo736NHD58+fXyzG0qxfv95HjRrl\nXbp08Xr16nlWVpa3atXKL7jggiLf45lnnun16tWL20Zpx+JZtmyZZ2Zmetu2bRM+J92PLpuXknlG\nM7NNhDWC3gKeBZ5z9wUpz57SxMw6AnPmzJlDx44d0xrLihWQlwfR3bidO8MDD8ABB6QvLpHtwdy5\nc+nUqRNV4WddJB1+/PFHdt99d6677jquvvrqhM4p6+cm/zjQyd3npjbi8o1Z2RU4A/gKuBz42MwW\nmdldZtbLzKruYhLVTG4uzJgBV15ZONvt7NmhW+gvfwndRCIiIsl48MEH2bp1a7V4CihfwsmKu69y\n98nufgbQGDgUeAToCjwH/Ghmz5jZOWa2Z2ltSdlq1AhrB735JkSmNmDzZrjhBujQAV54IXQUiYiI\nJGLWrFncfffdjBw5khNOOKHcj4enU1JDgCNdVO+5+9/c/SCgCTAUWA/cDCwxs3lmdqOZtS61MSnV\noYfChx+GOyqZkeHQCxbA0UdDjx6Q5PIaIiKygxk+fDiXXnopHTt2TOipsKokJc8rufv37v6gu58M\n5AK/A2YAxwKnpOIaO7JatWD4cJg7Fw4+uLB81izo2hVOPDEkMCIiIiWZNWsW69ev5+WXX67aixbG\nkfKHq919s7vPcvfL3H1/4KZUX2NH1b49vPMOTJ4M0TM8P/kk7L8/HHccvPaauodERGT7ktA8K2Y2\nMIm23d0fcfdNSZwrJcjIgFNPhX79wtNB118Py5aFBGXatLB16AAXXwynnKJZcEVEpPpLdFK48Um0\n7YQBuFIJatYME8j94Q9hMcQxY+Drr8OxDz8M5cOGQf/+YVXnDh0KnywSERGpThLtBmqRxNYy1cFK\ncTk5cNllsGgRTJwYHm/Ot2JFSGI6dYIDDwzrEH2pRRBERKSaSShZcfelyWyVHbwUqlkzTCT3wQfw\n+uuhCyh6Da558+DSS6FZs5C83HADfPKJxreIiEjVl/DaQCUxs/2A/CU4l7r7/Iq2KckzgyOOCNtP\nP4W1hR56KAzMzTd3btj+8hfYe2846qiwde8ODRumL3aRbWmBHqETSVi6f16STlbM7DhgFNA8pnwx\ncLG7T6tYaFJRDRvCOeeE7dNP4fHH4d//hjlzCut88UXY7r03DN7t1Am6dYPDDgtzvCS4HpdItZGb\nm0vt2rWr1eydIlVB7dq1yc3NTcu1E14bqMhJZr2BacBS4D4gP+VqCwwh3Gnp4+4zUhRnpatKawNV\nti+/hKefhqeeCjPkbirlea1WrULS0rkzdOwYxr6kaKFUkbT58ssvWbFiRbrDEKlWcnNzS5z1trLX\nBko2WXkHqAUc4e5rY47lAG8C6939kJREuQ3sSMlKtF9+gTfeCIsmvvRSGNtSmowMaN06PF3Url2Y\n32X//aFFi3BMRER2PJWdrCTbDXQAcHVsogLg7mvNbDwwMpmGzawzcCbQjdDF9CPwLnCtu3+ewPn1\ngVuB44HawPvAJe7+YTLxbO/q1IFjjgkbwPLl8Pbb8NZb4XX2bNi4sbD+1q1httzY7svatUMS06pV\n2PbdN7y2bAmNGumxaRERSV6yycp6YOdSju8cqZOMKwiLJD4O/JewaOJQYK6ZdSltAK+ZGWFRxfbA\nLYRE53zgVTPr6O5fJBnTDqNRIzj++LABrF8PH30U5m7JH5g7b17xrqN160KdD+OkhNnZ0Lx5uPvS\nrBnsuWfYmjYNr02ahDoiIiLxJJuszAT+bGYz3P2d6ANm1gW4EHgxybZvB/LcfXNUm48B84ArgdJm\n0z0ZOAQ40d2fipz7OPAZcD2gEXXllJUV1h/q2rWwbMOGMGD3k0+KbosWwZYtxdv49df4d2Oi1a8P\nu+9euDVqVHTbdVfIzYVddgl11eUkIrLjSDZZuRx4B3jTzN4HPo2UtwYOBpYT7pCUm7u/G6dsoZl9\nQhjAW5oTgWX5iUrk3BWRZOd0M6up6f8rrlYtOOCAsEXbtAkWL4bPP4fPPguvS5aEsiVLwl2akvz8\nc9j+97+yr1+jBuy8c9gaNiy61a8PDRqE1/ytXj2oW7fwtW7dMC+NiIhUD0klK+6+2MwOAK4CjgFO\njRxaCtwJ3OTuy1MTYoHdgI/LqNMBiDew533gbGBf4JMUxyURNWuGsSr77gu//33RY+7w/fewdGlY\nFuDrr+Grr8L23XeF29pio6CK27IFfvghbMnaaacwXid/y8kJ425ycgr387fs7MLX/C0rq3C/Vq3w\nPisr7MfbatRIPlYRkR1d0vOsRJKRYZGtUpnZAGAP4Noyqu4OvBan/LvIaxOUrKSFGTRuHLYuXUqu\nt2ZNSFp++CEM9o3efvyxcFuxAlauDPWTsXFjOH/lyuTOL68aNUKCFG+rWbPwtbQtM7P4a7ytRo3i\n+zVqFG6x72O3jIz476PLMzJK3mKPm5VeP/p4/r4GZItItArPYFvZzKwNcDfwFvBwGdWzgQ1xytcD\nFjkuVVh+N82++yZWf/NmWLUqzNb700+F3UnR25o1sHp1eM3fX7s2PLb9yy+hbN26yl16YMuWMHbn\n118r7xrbG7OiyUvs+9jXeFv0segkqKT6JW2lnZPosbL2Sysr6TW2LF690tqoyLFErhfvfSJtlfS+\nImXb+rx4UnlesvUq63oVudOdiIrMYNsMOIOwYGFDQjIQzd39uArEhpntBjwL/ASc7GVPCvMrYf6X\nWFmEVaD1q2I7k5kZBt5WdFJF9zCmZu3awi0/uVi3rvB1/fqwH/26YUN4zS/buDGURW8bNxZu+e83\nbSpavnlz2XHuSNzDtnVruiMRkXRLKlkxszzgocj5q4Cf41Sr0N+pZlYPmAHUAw5392UJnPYdoSso\nVn7Zt2U1MGzYMOrXr1+kLC8vj7y8vAQuL9WVWeEYlDTNJo17SFjyk5hNmwrfR+9v3ly4bdoU7tpE\nl23eXLQsf3/LluL7sdvWrSWXRb9Gb1u2FCYV0fXyy7ZuDfuxZdHvY+vGK4tOXmL3Y+vEq19anfw/\ng8pTXp5j+fsi24dJkS1avDQgdZKdwXYJsBo4yd0/S3lQZrWAlwgDZnu4+/sJnvcYIbFpElN+H5AH\n7FzS00A76gy2IrLtxUtiSior6TVeElRS/UTLKtJW7LF475M9rzLaiqe8MZSn/VSel2y9yrzeggVz\nycurejPY5gK3VFKikgE8BnQBji0pUTGzxkB9YKG758/uMRU40cz6ufuTkXq5wEnAND22LCJVQfR4\nE5HtQbw5tlIp2WTlPSD+akYVNwroS1goMdfMTo8+6O4TIrs3ESaIaw58GSmbClwEPGhm+wMrCDPY\nZgDXVVK8IiIiUomSTVYuAp43s9nuPjWVAQEHEsa79I1ssfKTFQeKDL1z961mdgxhbaChhKd/3gcG\nJrKukIiIiFQ9yU4KN8/MrgEmm9la4Gsg9iaQu/uBSbTdPcF6g4BBccp/BoZENhEREanmkn0a6Hzg\nLsL8JV9Q2cOARUREZIeVbDfQ1cDbQJ/InQwRERGRSpHs2rX1gQlKVERERKSyJZusvAa0T2UgIiIi\nIvEkm6ycB/zWzC43s11SGZCIiIhItGSTlflAC+BGYLmZrTWz1TGbuohERESkwpIdYPsEFVz7R0RE\nRCQRyc6zcmaK4xARERGJK9luIBEREZFtIqFkxczyzMq/7JYFeeUPS0RERCRI9M7KHcBnkad/WpRV\n2cz2MbOrgYXA6IoEKCIiIju2RMestCQsXngJcKOZLQHmAouBnwADGhKeEOoM7An8CIxByYqIiIhU\nQELJiruvBf5uZjcTVkI+DjgU6EdIVCA8HfQFYcK4p4Fn3H1TyiMWERGRHUq5ngZy983AU5ENM6sB\n7Bw5vNLdY1deFhEREamQZOdZASCSnPyQolhEREREitGjyyIiIlKlKVkRERGRKk3JioiIiFRpSlZE\nRESkSlOyIiIiIlWakhURERGp0pJ+dNnMegFnEWa3bUjh5HD53N33rkBsIiIiIsklK2Z2GXAT8D3w\nPjAvlUGJiIiI5Ev2zsqfgZlAb02pLyIiIpUp2TErDYGpSlRERESksiWbrLwPtE5lICIiIiLxJJus\nnA/0M7P+qQxGREREJFayY1amRM59xMzuAb4GYldcdnc/sCLBiYiIiCSbrKwEfgQ+T2EsIiIiIsUk\nlay4e7cUxyEiIiISl2awFRERkSqtIjPY1gAGAL8HmkWKlwLTgQnuHjuGRURERKTckrqzYmb1gbeA\nfwE9gZqR7SjgQeBNM6uXqiBFRERkx5VsN9DfgU7AUGBXd+/o7h2BRsAFQOdIHREREZEKSTZZOQEY\n6+5jo2exdfdN7n4PcA9wYioCFBERkR1bssnKLsCnpRz/H7Bzkm2LiIiIFEg2WVkIHFvK8WOBL5Js\nW0RERKRAssnKWKCnmT1nZj3NrHlk62VmzxIG2t6dujBFRERkR5XspHBjzawRcCXQK+bwJmB4ZOxK\nUswsB7gcODiyNQTOdPeHEzj3DMITScXCBnZ39+XJxiUiIiLbXtLzrLj7dWZ2N/A7is6z8rK7r6hg\nXLnAXyLtfQR0K294kfOXxJSvqmBcIiIiso0lnawARJKSySmKJdq3QGN3X25mnYAPkmhjhrvPTXFc\nIiIiso0llKyY2V4A7v5l9Puy5Ncvr8jj0BXurjGzOsA6d99a0bZEREQkPRK9s7IEcDPLdveN+e8T\nOK9GknFVlAGvAnWAjWb2AnCJuy9MUzwiIiKSpESTlcGE5GRTzPuqaB1hgO0sYDVhpt1LgLfMrKO7\nf5PO4ERERKR8EkpW3H18ae+rEnd/HHg8qmiamb0IvA5cA5yflsBEREQkKckuZPgvM+tSyvGDzexf\nyYeVWu7+FvAe4cklERERqUaSfRroTOBlQgIQTwvgDEJ3UVXxFbBvWZWGDRtG/fr1i5Tl5eWRl5dX\nWXGJiIhUG5MmTWLSpElFyn7++edKvWaFHl0uRRPg10pqO1ktgR/KqjR69Gg6duy4DcIRERGpfuL9\nAT937lw6depUaddMOFkxs+OA46KKhphZvG6VBoTulmTmRikXM2sM1AcWuvuWSFlu7KR0ZtabMND2\njsqOSURERFKrPHdW9gNOjuw70IWQAERzYC1hMOvFFQnMzP5ESHz2iBQda2Z7RvbHuPsa4CZgINAc\nyJ/T5W0z+xCYDfwciXEQYTbcGysSk4iIiGx7CScr7n4jkV/2ZrYVOMvdJ1ZWYMClQP7kcw6cENkA\nHgHWRMpjJ3ybDPyesJhibeA74J+E9YrK7AYSERGRqiXZhQyTXa25PNdokUCdQYS7JtFlfwX+Wllx\niYiIyLZV6UmHiIiISEUknayY2TFm9pKZ/Whmm81sS+yWykBFRERkx5TspHAnAtOB3QhjRDKASZH9\nX4H/AsNTFKOIiIjswJK9s3IV8D7QAfhbpOxf7n460A7YHVhc8fBERERkR5dssrIfMDkyt8nmSFlN\nAHdfAowFrqhwdCIiIrLDSzZZWQdsBHD3VcAGwt2UfN8TptwXERERqZBkk5VPCXdX8n0E/MHMMs0s\nC+hP4SRtIiIiIklLNll5CjjOzGpF3v8d6AasIqy/cwRhdlkRERGRCkl2UrjbgNui3k83s25AP2AL\n8Ky7z0pJhCIiIrJDK3eyErmb0gtY4u7/zS939zeAN1IYm4iIiEhS3UAbgceBQ1Mci4iIiEgx5U5W\n3N2Bz4Hc1IcjIiIiUlSyA2xHAheYWetUBiMiIiISK6kBtkBX4EfgYzN7FVhCmGY/mrv7n5MPTURE\nRCT5ZOWCqP0eJdRxQMmKiIiIVEiyjy4nvVqziIiISHkku+ryXmaWXcrxbDPbK/mwRERERIJk75As\nBk4o5fixaNVlERERSYFkkxUr43hNYGuSbYuIiIgUSHjMipnVAxpEFe1SQldPA+A04LsKxiYiIiJS\nrgG2w4C/RvYduCOyxWPAtRWIS0RERAQoX7LyIvALIRG5BZgEzI2p48BaYI67z05JhCIiIrJDSzhZ\ncfd3gHcAzCwHeMLdP66swEREREQg+XlWrk91ICIiIiLxJDuDLWbWEMgDWgINKf6EkLv7WRWITURE\nRCS5ZMXMegFTgRxgNfBTnGpegbhEREREgOTvrNwOLAP6ufu8FMYjIiIiUkSyk8LtA4xRoiIiIiKV\nLdlk5XOgbioDEREREYkn2WTlWuB8M2ueulBEREREikt2zEoP4AdggZm9BHwFbImp4+7+54oEJyIi\nIt6iCCYAACAASURBVJJssnJB1H6fEuo4oGRFREREKiTZSeGS7T4SERERKRclHSIiIlKlJT2DLYCZ\ndQW6A42Ase7+uZnVBtoAn7n7LymIUURERHZgSd1ZMbOdzOxJ4C3g78CFwJ6Rw1sJKzRrvIqIiIhU\nWLLdQCMIA2vPA1oTtS6Qu68HHgeOq3B0IiIissNLNlnJA+5x9/uAlXGOLyAscJgUM8sxs+vN7Hkz\n+9HMtprZwHKcX9/M7jOz5Wb2i5nNNLMOycYjIiIi6ZNsstIIKG2q/S1A7STbBsgF/kIY+/IR5VgU\n0cwMeA44DRgDXMb/t3ffcXLV9f7HX5/t2SRseiMJCb0FUjW51IsU0Ws0NMmlg6CgXAUE5AcBRLlI\nE9QriHRFkSIQEKLUCAQhDQyEkkISIIX0sMlutn5+f3xn22R2d3Zmdmd29v18PM7jnPme7znnMyf7\n2P3knG+B/sBMM9stiZhEREQkDRJtYPspIZFozkHAkgTPDbAKGOTua81sHDCnDceeCEwCjnf3JwHM\n7DFgEfBT4NQk4hIREZEOluiTlT8D3zWzSY3KHMDMzgVOAv6QaFDuXuXuaxM8/HhgTV2iEjnfeuBR\n4Jtmlp9oXCIiItLxEk1WrgfeAF4FXiEkKreZ2SfAXcDfgdtSEmHbjQHmxyifTXg1tWfHhiMiIiLJ\nSChZcfdK4KvAWcDHwIdAIbAAOBP4hrtHzxXUUQYDq2OU15UN6cBYREREJEkJDwrn7g48FFkySTeg\nIkb5dkIX624tHl1a2g4hiYhA+LUJjjfZrtvXeDt6X2t1W6ofb1ky54reF32eZI6LVSdWzC0d29r1\nYh3Xnuduy7GpiqE9z7Vo/aK4rpeohJIVM+sDDHX3Bc3sHwV85u6bkgkuQeWEpzzRigivq8pbPPrd\nd+Gww9ohLJGWuTtVtVVU1VRRWVNJdW11/eeq2iqqa6tDWU3DdnNLjdc0bNfW1JdFb0eva722SVmt\n1+6wXeu1DQu1TT7vsD+yOB6zvH6/+w5168oaf677I99cvcZ1ouu2Vg+Iua/u3yaZfSJZb1X7nj7R\nJyu3EQaDm9jM/rsIY62ck+D5k7Ga8CooWl1Zi7f0oltuoeT555uUTZ06lalTp6YmOslI7k55dTnb\nKrexrWobZVVllFWVUV5VHtbVYb29ejvlVeVhXR3WFdUVbK/eHpaa8LmipqJ+XVlTSUV1WFfWVNaX\n1SUldUtN2t6cioi0wbvsOHjJ9va9ZKLJyhHAnS3sfwb4XoLnTtY7wMExyicCZYQuzM26bZddGPv0\n0+0Rl7SD6tpqNpVvYtP2TWzevpnN2zezZfsWtlRsqV+XVpTyRcUXlFaWhqWilK2VW5ss26q2pfur\nSJQcy8EwzKzJthE+51hOzH2N17GOA2LWbWnd0jHx7mttG6g/tm47el9rdWPVb+kcLV2nteOaOz76\nHK1dL9a5Yh3XXJ32Ple814vr3DHKYh6b4PkTOU9C5xq34/51S9bxxEVPxHW9RCSarPQH1rewfwNh\n4Lh2ZWaDgBJgSaMGvY8Dx5vZce7+RKReP+AE4Gl3r2rxpO+/344RS0vcndLKUtZsXcPabWubLOu2\nrWN9+Xo2lG1gQ/mG+vXWys4xV2ZeTh4FuQX1S35OPoV5hfXb+bn5TbZjrfNy8sjLyasvy8vJI9dy\n67frllzLbdjOya0vy83Jrd9Xtx1rnWM5cW3XJwxYk7Jcy8XMyLXcHRKKxvXqPkcnH43/CItI5zB/\n/nyeIPOSldWELsLNGQesS/DcAJjZ94FewM6RoslmVjdZ4q/dvRT4BXA6MAL4JLLvceBHwP1mth8h\nqbqA0PPp2lYvvH49rFoFQ9RpKJXcnTVb17Biywo+++IzPvviMz7d8imflX7GqtJVrC5dzeqtqymr\nKuuQeApzC+lR0IMeBT3oXtCdHgU9KM4vpnt+d7oXdKc4v5jivGKK84vplt8trPO60S2/G93yulGU\nV1S/XZhXSFFeEUV5RRTmFlKYV9hkXZBbQG5Obod8LxGRbJRosvIU8H0zm+HuTd6ZmNk3CV2aW3pN\nFI8fA8Mj2w5MiSwAfwRKI+W1jQ9y91ozOxa4GbiQ0PtnNnC6uy+O68pz5sA3NQ9jW1XWVLJs0zIW\nb1zMog2LWLJxCcs2L2PZpmWs2LKC7dWpeamZa7n0Le5Ln2596F3Um97deod1UW9KikroVdSLksIS\nSopKKCksYafCnehZ2DOsC3rSo6AH+bkaG1BEpLNINFm5FjgSeNLM/g28FynfHziQ0Lj2mmQCc/eR\ncdQ5i5AYRZdvAc6LLG03d66SlRZUVFfw4foPWbhuIQvXLgzrdQtZtmlZwo1EexX1YnCPwQzuOZhB\nPQYxsPtABnQfUL/0L+5Pv+J+9C3uS0lhiV4TiIh0IQklK+6+xcwmApcBxxHagwAsBX4G3OzunbfF\n4py2TEWU3bZXb+edNe8wf/X8+uW9te9RVdty05/GivOLGdFrBCN7jWSXkl0YVjKMYTsNY+hOQxlW\nMozBPQbTLb/l4W9ERKTrSmZQuG2EpydJPUHJSHPmgDt0wf+9r922llmfzGLWp7N449M3mLd6HpU1\nla0eV5xfzF5992LPvnuyR589wrrvHuzae1f6F/fXkxAREUlYwslKVtu4EZYvh5Gtvonq9LZWbuXV\nFa/y4scv8sLHL/De2vdarJ9jOezdb29GDxrNqAGj2K//fuw3YD9G9BpR3y1UREQklRJOVsxsH0J7\nkV2B3rBDR213968kEVt6zZmTtcnKis0rmP7RdJ768Cle/+T1Fl/p7Nl3TyYNncSEIRMYO3gsBww8\ngO4F3TswWhER6eoSHW7/NOB+oAr4CIg1rH7nfu4/dy6cdFK6o0iZj9Z/xKMLH+Wpj55i/upYk1KH\npybjh4zn8F0O56DhBzFp6CT6d+/fwZGKiIg0lUxvoLeBY929pcHhOq8saGS7qXwTjyx8hAf//SBv\nfvZmzDq79d6No3c7mqN2PYrDRxxO7269OzhKERGRliWarAwBbsnKRGXAAFi7FubNg9payOlc7TDc\nndc+eY3fzvkt0z+cTkXNjhNQjxs8jil7T2HKPlPYp98+avwqIiIZLdFkZQEhYck+++4bkpXSUli0\nCPbeO90RxaWqporH3n+MX/7rl8xbPW+H/aMGjOKMA8/gxP1OZHjJ8BhnEBERyUyJJisXA49FRrB9\nI5UBpd2++8LMmWF7zpyMT1a2VW7jjjl38Ku3fsXK0pVN9vUr7scpo07hjAPPYPSg0XqCIiIinVKi\nycrlwBbgNTN7nzAvT/TQpe7unW8Y2P32a9ieOxdOOy19sbSgqqaKe+bfw3WvXsearWua7Bs7eCwX\nT7yYE/c7kYLcgjRFKCIikhqJJisHEObl+QToAewbo44nGlRa7bNPw3YGNrKt9VoeXfgoV718FUs3\nLa0vN4zJe03m4kkXc8jwQ/QURUREskaiw+2PSHEcmaOkBHbbDZYuhbffhupqyMuMsfMWfL6Ac585\nl9krZzcpP36f4/n5ET9n736Z/cpKREQkEZ2rq0tHGT8+rLdvh4UL0xsLYeLAaS9PY9zvxzVJVP5z\nxH/y1nfe4vGTHleiIiIiWSupRwZmdhjwdWCXSNEK4Fl3/2eygaXVhAnwyCNhe84cOPDAtIXyxqdv\n8J2nv8MH6z+oL9un3z7c/tXbOWrXo/S6R0REsl6iI9gWAA8D3yKMVLs5sqsXcImZPQlMdff4p+bN\nJBMmNGzPnQvf+U6Hh1BdW820l6dx46wb8Ujzn7ycPK44+AquPORKCvMKOzwmERGRdEj0yco1wBTg\nFuBWd/8cwMwGAJcAlwJXA9NSEWSHGzMmzLjsnpZGtuvL1jP1r1N58eMX68smDJnAPZPv4YCBB3R4\nPCIiIumUaJuV/wYedPfL6hIVAHdf6+6XA38AMrPPbzx69mzoFbRgAZSXd9il562ax7jfj6tPVPJy\n8rjpyJv41zn/UqIiIiJdUqLJymDgrRb2vwUMSvDcmWHSpLCurobXXuuQSz7wzgMcdN9BfLLlEwAG\ndh/Iy6e/zKUHXUpuTm6HxCAiIpJpEk1WPgMOb2H/YZE6ndcxxzRsz5jRrpdyd6555RrOmn5W/Vw+\nk4ZOYv5353PILoe067VFREQyXaLJyoPASWb2OzPby8xyzSwnsn0ncCLwQMqiTIejjoLcyNOMdkxW\n3J0rXrqC6169rr7sgvEXMPPMmQzpmZ3TL4mIiLRFog1s/xfYDTgPOBeojZTnEHoHPRip03n16hVe\nBb3+Onz0ESxbBiNHpvQS7s7F/7iY29+6vb7s9mNu54cTf5jS64iIiHRmCT1Zcfcadz8TGA1cCdwT\nWa4ERrv7We5e28IpOodjj23YTvHTlVqv5QfP/aBJonLn1+9UoiIiIhIlqUHh3H0BsCBFsWSeY4+F\nK68M2zNmwAUXpOS07s73/vY97p5/NxDm9bln8j2cPebslJxfREQkm8T9ZMXMiiJtVC5spd7/mNmd\nkYHjOrfRo2FQpFPTyy+H4fdT4IbXb6hPVHIshz9M+YMSFRERkWa05TXQecCZwLOt1HsWOCtSv3Mz\ng69+NWyXlaWkC/P0D6dz5cvhaY1hPDTlIU494NSkzysiIpKt2pKsnAT81d0/bqmSuy8FHgemJhNY\nxkhhu5V3P3+XU59sSEx+fsTPmToqO26TiIhIe2lLsjIKeD3OurOA/dseTgY66ijIidymJJKV9WXr\nmfyXyWyt3ArAyfufzBUHX5GKCEVERLJaW5KVAqAyzrqVkfqdX+/eDaPZfvghLF/e5lNU1VRxwqMn\nsHxzOHbc4HHcO/lezZgsIiISh7YkK6uI/2nJ/pH62SHJV0GXPH8J/1zxTwAG9RjEUyc/RXF+caqi\nExERyWptSVZeBE6PzKzcrMj+04EXkgksoySRrLyy7BV+M/s3ABTmFvLkt59k6E5DUxmdiIhIVmtL\nsnIjUAS8bGZfjlUhUv5SpN7NyYeXIUaPhoEDw/ZLL0FFRVyHlVWVce4z59Z/vuXoW5g4dGJ7RCgi\nIpK14k5WIr2ATgKGA2+Y2WIze8LMHoysFwFvACOAkyO9grJDTk7TLsyvvhrXYVe/cjVLN4XbcMjw\nQ7hgQmoGlRMREelK2jTcvrs/CxwA/J7w9ORbwGmRdTFwN3Cguz+T4jjTr42vgmavnM1tb94GhNc/\n90y+hxxLdN5IERGRrqvNfz3dfbm7n+/uw4ASYBhQ4u5D3f17rY3D0mkdfXRDF+YnnoDa5qc+qqyp\n5Jynz6E2Mj3STw//KXv23bMjohQREck6Sf1X391L3X2lu5emKqCM1bs3HHNM2F6xIgy/34wbXruB\n99a+B8DYwWO55D8u6YgIRUREspLeS7TF2Y3m77nvvphV3lv7Hte/dj0AuZbLvZPvJS8nqfkiRURE\nujQlK23xjW9A375h+4knYNOmHapc9sJlVNVWAXD5QZczetDojoxQREQk6yhZaYvCQjjttLBdUQF/\n/nOT3W9+9iYzloTGt8NLhjPtsGkdHaGIiEjWychkxcwKzOxGM1tpZmVm9qaZHRnHcWeYWW2Mpaa1\nwezi1sKroGtnXlu/fdUhV1GUV5SSS4qIiHRlmdqY4kHgOOA2YAlwJvCcmR3u7m+0cqwD04DlUeWb\nUxLZqFEwfjzMnQvz58M778Do0fzr03/xj6X/AGBErxGcOfrMlFxORESkq8u4ZMXMvgR8G7jE3W+L\nlP0ReA+4CTg4jtP83d3nt1uQ55wTkhUIT1d+/WuumXlN/e6rDrmK/Nz8dru8iIhIV5KJr4FOAKoJ\nA8wB4O4VwL3AJDPbOZ6TmFkPs3Yahe3kk6Eo8ornoYeYteRlXvg4TIW0a+9dOf3A09vlsiIiIl1R\nJiYro4FF7r41qnx2o/0tMWAm8AVQZmbTzWz3lEbYqxccf3zY3rSJa578Yf0uPVURERFJrUxMVgYD\nq2OUryYkIkNaOLYMuB+4gDAFwI3AV4BZ8T6RiVukoe1rw+GlsjAA3G69d+O0A09L6WVERES6uoxr\nswJ0A2JNa7y90f6Y3P0x4LFGRU+b2fPAq8CVhCQmNQ4/HEaO5JpDl9UXTTt0mgaAExERSbFMfLJS\nDhTGKC9qtD9u7j4LeAtotetzm+Tk8K8zv8IrI8PHPWp7c8oBp6T0EiIiIpKZT1ZWE/tVz+DIelUC\n5/wUiGsmwYsuuoiSkpImZVOnTmXq1Kk71L17j1JYFLav+Mc28n6wCfr3TyA8ERGRzuHhhx/m4Ycf\nblK2ZcuWdr2muXu7XqCtzOwm4EdAn8aNbM3s/wE/A4a7+8o2nnMO0MPd92mhzlhg3rx58xg7dmyr\n5yytKGXwrYPZVrWNnbbD6luh+EeXwk03tSU0ERGRTm/+/PmMGzcOYFx7DB2Sia+BHic88TmvrsDM\nCggDw71Zl6iY2SAz28vMchvV6xd9MjP7GjAOmJHKIB9d+CjbqrYBMPX9HIqrgN/+FtauTeVlRERE\nuryMS1bcfTahkewNkSH3zwVeAXYBLmtU9RfAB0DjXj5vmNkjZnapmZ1nZncBTwErgBtSGed97zQM\ntX/2bieEjbIyPVkRERFJsYxLViJOA24HTgV+BeQCX480lq3jQG3UcX8BdgeuAH4NHA3cBXzJ3del\nKrgP13/IG5+GUf/3H7A/E358W8MgcXfcAWvWpOpSIiIiXV5GJivuXunul7v7zu5e7O4T3f3FqDpn\nuXueu3/SqOxqdx/n7n3cvcjdR7r7halMVADuf/v++u2zR5+NDRkC3/teKCgv19MVERGRFMrIZCWT\nVdVU8eC/HwQgPyefUw84Ney4/HLoFhkC5s47YXWsce1ERESkrZSstNGMJTP4fNvnAEzeazL9u0e6\nKg8aBOefH7a3b4cbb0xThCIiItlFyUob3fd2o4a1Y85uuvOyyxqervzud/Dxxx0YmYiISHZSstIG\na7au4W+L/gbAzj135pjdjmlaYeBAuPDCsF1RAeedBxk2jo2IiEhno2SlDR5a8BA1XgPAGQeeQW5O\n7o6VrroKhg8P2y+9BPffv2MdERERiZuSlTi5e5NXQGeNOSt2xZ494a67Gj5ffDGsSmSGABEREQEl\nK3FbvHExH6z/AICDhx/M7n12b77yV78Kp58etrdsgQsu0OsgERGRBClZidOMxQ2j9U/ec3LrB/zy\nlzBgQNiePh0ef7ydIhMREcluSlbiNGNJQ7Jy7B7Htn5A377wf//X8PkHP4ANG9ohMhERkeymZCUO\nZVVlzFw+E4ChOw1lv/77xXfgCSfAlClhe+1a+P739TpIRESkjZSsxGHm8plU1FQAcOzux2Jm8R1o\nFmZiLikJnx95BG6+uZ2iFBERyU5KVuLQuL3KsbvH8QqoscGD4b6GXkT85Cfw7LMpikxERCT7KVmJ\nQ117lbycPL6y61fafoLjjoNrrw3b7jB1Krz/fuoCFBERyWJKVlqxeMNilm5aCoQuyzsV7pTYiaZN\nC21YAEpLYfJk2LgxRVGKiIhkLyUrrWjSC6itr4Aay8mBBx6A0aPD56VL4aSToKoquQBFRESynJKV\nVqQsWQHo3j2MudI/MlPzSy+FweOqq5M7r4iISBZTstKC8qry+i7LO/fcmf0H7J/8SYcPhyefhPz8\n8Pkvf4GTT4bKyuTPLSIikoWUrLRg5vKZbK/eDrSxy3JrDjoInngCCgrC57/+NbRnqahIzflFRESy\niJKVFrR51Nq2+K//gqefhqKi8PmZZ+Bb34Ly8tReR0REpJNTstKCxl2Wj9z1yNRf4JhjwpgrxcXh\n89//HpKYTZtSfy0REZFOSslKM5ZsXMKSjUsAOGjYQYl3WW7NEUeEJKVHj/D55ZdhwgR49932uZ6I\niEgno2SlGUmNWttWhxwCzz8fJj+E0K154sTQ+FZERKSLU7LSjBeXvVi/nfL2KrFMmgTz5sHYseFz\nWVkY6fbiizUWi4iIdGlKVmJwd2avnA1Ar6JejBowqmMuvMsu8PrrcOaZDWW33QaHHqrh+UVEpMtS\nshLDqtJVrNm6BoDxQ8anrstyPLp1CxMf3nFHw1gsb74JY8bAz36m8VhERKTLUbISw5xVc+q3xw8e\n3/EBmMH558Orr8Iee4Syykq4+moYNw5mz+74mERERNJEyUoMc1fNrd+esPOE9AUycSL8+9/wk59A\nbm4oe++90L7lO9+BTz9NX2wiIiIdRMlKDI2frEwYksZkBcJroRtugDlzwqsggNpauPfe8NTl0kth\nw4b0xigiItKOlKxEcff6JysDug9g6E5D0xxRxJgx8NZbcNNNUFISyioq4JZbYNdd4ec/h40b0xuj\niIhIO1CyEmVl6Uo2loc/+hOGTOjYxrWtyc8PT1I+/jisCwtD+RdfwLRpMGwYfP/7sHhxeuMUERFJ\nISUrUd5f29BFePyQNDSujUefPuEJy5IlcM45kBP5ZywrC72I9toLvvlNeOEFqKlJb6wiIiJJUrIS\n5f11DclK2turtGboULjnHvjoI7jwQujePZS7h0kSjz46jN1yxRXw4YfpjVVERCRBSlaiLFy3sH47\nY5+sRNt9d/j1r0PvoBtvhJ13bti3ciX84hewzz6hd9Gtt4bXSCIiIp2EkpUoH6z/AIBhOw1jYI+B\naY6mjXr3hssug2XL4LHHwgzOdV2eITTQ/fGPYbfd4IAD4JprwhD/tbXpi1lERKQVSlailFeVA2ke\nXyVZ+flwwgnwzDPhycovfxmSk8befReuuw7Gj4eBA+Hkk8MrpRUr0hOziIhIM5SsNCMtI9e2h4ED\n4aKLwuByH3wQxmz58peb1lm/Hh55BM49F0aMgJEj4dRTQ2Pdd95RI10REUmrvHQHkKk69ZOV5uy9\ndxgN9yc/gVWrYPp0eP55ePnl0P25zvLlYfnTn8LnHj3COC/jxoVZoceODT2O8vTjIyIi7c/cPd0x\nZAQzGwvM4zxgCGy8bCO9u/VOd1gdo7oa5s4NXZ1feinMPVRe3vIxRUWh0e5++4Vl//1DMjRihJIY\nEZEuZv78+YwbNw5gnLvPT/X5M/KvipkVAD8DTgV6AwuAq9z9xTiOLQFuBr4FFAOzgUvc/e14r797\nn927TqICIbmYODEs06ZBVVV4/fPGGzBrVpj1OXoeou3b4e23wxJ9rl13DVMB7LFH2B45MiQxI0aE\npzQiIiJtkJHJCvAgcBxwG7AEOBN4zswOd/c3mjvIwnCzzwGjgJuADcAFwEwzG+vuS+O5eKfpstxe\n8vNhwoSw/PCHoWzdupCYzJ8fehAtWBAGpYvuSVRdDYsWhSWWvn3DSLvDhoVxYoYNgyFDYPDghnWf\nPmHmaRERETIwWTGzLwHfJjwNuS1S9kfgPUICcnALh58ITAKOd/cnI8c+BiwCfkp4UtOqjB8MLh36\n9w+DzB19dEPZ9u1hQLqFC8NSl6QsXtz8a6QNG8LyzjsAPAxMja5TUBCuN2BAw9K/f0h0+vUL6759\nQ1LTu3dYFxcrwYnTww8/zNSpO9x1aUe65x1P9zy7ZFyyApwAVAN31xW4e4WZ3Qtcb2Y7u/vKZo49\nHlhTl6hEjl1vZo8Cp5hZvrtXtRaAkpU4FRXBgQeGpbHa2tCAd/Hi0FB32bKG9YoVYV+kh1HMZKWy\nMnS5XtncP3MM+fkhcSkpgV69wrpu2Wkn6NmzYd2zZ3gd1Xjp3j0kPN27h5muszjx0S/xjqd73vF0\nz7NLJiYro4FF7r41qnx2o/3N/RUbA8Rq2DMbOBfYE1gYY389M2PM4DHxRys7yskJr3iGNjNjdU0N\nrFkDn30WJl484wxYvbrpsm4drF0bf7fpqqpQf+3a1HyH4uKGpVu3hnW3biFJi14XFYWJJevWdUtB\nwY7bBQVNl/z8hnWsJTc3q5MnEZHWZGKyMhhYHaN8NWDAkFaO/WczxxI5tsVkZWSvkfQoUCPQdpWb\nG6YE2Hnn0E7lwgtj16uthc2bG5KQuldIGzaEsWE2boRNm5ouW7aEbtjJ9nIrKwtLpsjPD42Xo9ex\nltzcHbdzcxuWefPguOOaljVecnJif87JaXm7uTKz+D/Xbcezjt6O9Tnefa0t0LZ9dWV16/Ly8FQx\n1jHRdVvbjq7f0rlaqt+WfbE+J3pcS2UizcjEZKUbUBGjfHuj/Ykca60cC8C+/fdtrYp0lJyc0B6l\nT5/QLTpetbWwdWtIXLZsgdLSkMCUljZsb9sW6tQt27btuJSXh6WsLKyrWn2D2H6qqsLSWpfyeD35\nZOt1JLVGjEh3BJkv0SQnVllNTXhimey5Ek20kknQ2vP8iSaJrR1XXZ3YeeOUiclKOVAYo7yo0f5E\njvVWji0C6LutL/Pnp7yLuDRjy5Yt7X+/CwtDA93+/ZM7T3V1aE9TUdGwVFY2lNVt1y1VVQ3b1dUN\nCUfd58ZlVVUNZXXl1dXhF25dWfR23efa2qZldUszthD7Xam0H93zOEU/EU3iCekWYH46/4PRxXzQ\nsFnUfK3EZdygcGb2PDDE3fePKj8CeBH4hrs/28yxiwjtXf4rqvxsQoPdA9w95msgM/tv4E8p+Aoi\nIiJd1Snu/udUnzQTn6y8AxxuZj2iGtlOJDwdeaeVY2N1bZ4IlBG6MDfnH8ApwHIaXjmJiIhI64qA\nEYS/pSmXiU9WvgS8CfzY3X8ZKSsgjLOyzt0PipQNAkqAJe5eEyk7idAb9kR3fyJS1o+QpMxw91M6\n+vuIiIhIcjIuWQEws0cIw+XfTsMItuOBI9x9VqTOA8DpwAh3/yRSlgO8DuwH3AKsJ4xgOwyY4O6L\nO/SLiIiISNIy8TUQwGnsODfQ1+sSlQgHmoz17u61ZnYsYW6gCwm9f2YDpytRERER6Zwy8smKiIiI\nSJ2cdAfQ3syswMxuNLOVZlZmZm+a2ZFxHltiZr83s7VmttXMXjYzDW/bikTvuZkNMrNfRO7zF2ZW\na2aHdkTMnV0S9/wIM7vXzD4ys21mttTM7o60CZMWJHHPDzGz6Wb2iZmVm9lqM5thZv/REXF3wfv/\n9wAACstJREFUZsn8Po86z92R3y9Pt0ec2SSJn/MzIvc4eqkxswFtjSPrkxXCDM4/Av4I/A9h3qHn\nWvvF0GgG55OBXwOXAv0JMzjv1q4Rd34J3XNgL8J9HkJ49afHfvFL9J7fCBwGPEF4dfowcBIwP5Ff\nKF1Movd8T6AGuJPQpu5mYCDwqpkd3dKBkvA9r2dm44EzaHncLWmQzD134CpCk4665TRgc5ujcPes\nXYAvEdq1XNSorBBYDLzeyrEnRY6d0qisH7AReCjd3y1TlyTveXegV2T7eMIv9EPT/Z0yfUnynh8c\no+yQyPmuS/d3y9QlmXvezPm6EaYFeS7d3y1Tl1Tdc2AWYdytZcDT6f5embwk+bvljMjv8LGpiCXb\nn6zEnMEZuBeYZGY7t3BszBmcgUeBb5pZfvuE3OklfM/dfZu7tz3jlmTu+esxyl4jJOX7pD7UrJHM\n75YduHs5sA7olcogs0zS99zMTif0Fr2yvYLMMin5OTezHpHeugnL9mQlnhmcm9PSDM7FhEe5sqNk\n7rkkJqX33My6Az0IXf8ltqTvuZn1NLO+ZraXmf0v4Y/oiymOM5skdc/NrAfwC+B6d0/R9OxZL9mf\ncwNmAl8AZZG2WrsnEkimdl1OlbTO4NxFJXPPJTGpvucXAfnAX5KMK5ul4p4/ChwT2a4E7gJ+npLo\nslOy9/wawkjmt6c4rmyWzD0vA+4HXiEkK+OAS4BZZjbW3Ve2JZBsT1bSOoNzF5XMPZfEpOyeR3pf\nXQ084u6xknUJUnHPLycMXjmM8H6/gJAkVqYiwCyU8D03sz0JjUO/7e6a3TB+Cd9zd38MeKxR0dOR\nuf9eJbyGu6AtgWR7spKuGZy7smTuuSQmJffczPYm9ApaAJybmtCyVtL33N0X1G2b2Z8Ir53vJzTu\nlx0lc89/RWgQ+lTKo8puKf197u6zzOwtoM3dzbM9WVlN7MdUgyPrVa0cOzhGeTzHdmXJ3HNJTNL3\n3MyGAc8DmwijRW9LXXhZKaU/5+5eFRnz43IzK4w0YpSmErrnZnYE4XXbFDPbpa6Y8PevW6Rso7uX\npjjebNAev88/JYE2n9newPYdYM9Iw6rG4p3BeWyM8nhmcO7Kkrnnkpik7rmZ9SEkKnnAMe7+ebtE\nmV3a4+e8mPBHtGeSsWWrRO/5sMj+JwndlZcBHxP+CH8lsn1WewScBdrj53xXQs+3Nsn2ZOVxwi/g\n8+oKLMzgfCbwZl0Dn8jIqXuZWW7UsQPN7LhGx/YjdOV6Wu89m5XMPZfEJHzPzawYmEH4n9LX3P3j\njgy8E0vmnvePPpmZ9SIMl/BJZIgE2VGi9/wlYAphctzGy3pgTmT7mQ76Dp1NMj/n/aJPZmZfIzS0\nndHWQLJ+biDTDM4dLtF7Him/ipCx70cYPfg+wv+EcPfrO+xLdDJJ/Jw/BUwmjJswM+q0W919egeE\n3yklcc/nAp8BbwFrgV0ixw4GTmo8tpM0lczvlhjnWga86+6T2znsTi2Jn/NFwNvAXGALIUk5C1gJ\nfMnd2/Z0Jd0j5LX3Qmhhf2PkBpUBbwJHRtW5nzDwzfCo8hLg94RfKKWEDH1Mur9Tpi9J3vNawqiH\n0Ut1ur9XJi+J3nNCIhjrftcAH6f7e2XyksQ9P58wLMLnhJ4WawivKP4j3d8p05dkfrfEONfHwPR0\nf6dMX5L4Ob8OmEcYYHJ75HfNb4D+icSR9U9WREREpHPL9jYrIiIi0skpWREREZGMpmRFREREMpqS\nFREREcloSlZEREQkoylZERERkYymZEVEREQympIVERERyWhKVkRERCSjKVkRERGRjKZkRaSLM7Nd\nzKzWzE5PdyyNWfCumV3RqOzaSKx90hlbIszs4cikcCLSRkpWRLJQ5A96a0uNmR0aOSQTJwn7b2Ao\n8H+NypzMjDUeNwLHm9modAci0tnkpTsAEWkXp0Z9PgM4MlJujco/cPd1ZtYNqOqo4OL0Y+Bhdy9N\ndyCp4O7vmNlc4BLgzDSHI9KpKFkRyULu/ufGn81sEmFa94ebqV/ZIYHFyczGAAcCF6U7lpaYWbG7\nl7XhkEeBa83sgjYeJ9Kl6TWQSBcXq82KmT1gZqVmNszM/hbZ/szMLojsH2VmL5nZVjNbbmZTY5y3\nxMxuN7NPzGy7mS02s8vMzKLrxvAtoAJ4rZn9vSMxbjKzzWZ2n5kVRV0/18ymmdmSyPWXmdn1ZlYQ\nVa/WzK6OEf9yM7uv0eczInUPNbM7zOxz4NPIvh6R77oscq3Pzex5MxsdddoXgB7AUXHcAxGJULIi\nIrE44ffDDGAFcCmwDPiNmZ0RKZ8DXAZ8ATxoZrvUHRx5rfQqod3JA8CFwOvADcCtcVx/EvCeu9fE\n2GeEJxTdgZ8AjxBec10TVe9e4KfAXOBHwEzgCiDm06UYmmsbcwewd+TcN0TK7gK+CzwGnA/cDJQB\n+0Qd+z5QDhwUZwwigl4DiUjzioA/uPtNEHqzAKsIScDJ7v54pPxF4ENCwnBd5NhLgJHAaHf/OFJ2\nt5mtBn5sZre6+8oWrr038GYL++e5+3l1H8ysH3AOIRnBzA4ATgd+7+7fi1T7nZmtAy4xs8Pc/Z9x\n3YUdrQe+4u6Nk5mvAXe7+2WNym6JPtDda8zsU2DfBK8t0iXpyYqItOTeug133wJ8BGyrS1Qi5YuA\nzcCujY47gfAKZ4uZ9a1bgJcI/0k6lJb1BTY1s88JTzIaew3oa2Y9Ip+/Fql3W1S9WwlPZr7eyvWb\n44SkJPqpy2bgy2Y2OI5zbAL6JXh9kS5JT1ZEpDnb3X1DVNkW4LMYdbcAvRt93gMYBayLUdeBAXFc\nv6W2LZ9Efa5LbHoDW4FdgFpgSZMLu39uZpsj+xO1PEbZZYTXXZ+a2TzgOcJTqWUx6hqdt/u1SFoo\nWRGR5sRqL9JSeePkIofQmPRGYicdi1q59gaaJj+JxADJJQW5zZSXRxe4+2Nm9iowBTia0O36cjOb\n4u7/iKrem9a/v4g0omRFRNrDUqCHu7+S4PEfEtq8JGoFIWHag/DqCgAzGwD0iuyvsylSRqN6+UA8\nr3TqufvnwO8IbWP6AW8DVwL1yYqZ5QLDgOltObdIV6c2KyLSHh4FJpnZ0dE7Il2am3tqUedfwP6R\npCERzxGesvwoqvwSwtOWZxuVLWXHNjTfpfknK02YWY6Z7dS4zN3XExojF0ZV35fQcHlWPOcWkUBP\nVkSkPdwMTAb+ZmYPAPMIXY0PAI4DRgAbWzh+OnAVcBjwYlsv7u4LzOxB4Dwz6w38E/gyoYfQE1E9\nge4hPA15nPDq6kDCq5xY7W1ivdLqCXwWOf7fhDYzRwHjgYuj6h4NbEvkO4l0ZUpWRLqOltpvxNrX\nXP3m6taXu3t5ZN6h/wecCJxGGI9lEXA1oUFu88G4zzezd4GTSPwP+zmEpyZnEgaZWwNcT0P36jp3\nE5Knc4BjCOPDHEXouRT9XWN99zLgt4REZArhifUS4Hx3/31U3ROAv7r7tkS+kEhXZTv2wBMRST8z\nO5UwieFwd/8i3fEkKzKa7VxgjLu/m+54RDoTJSsikpEiw/L/mzCZ4Q2t1c90kUH1cPcdpiYQkZYp\nWREREZGMpt5AIiIiktGUrIiIiEhGU7IiIiIiGU3JioiIiGQ0JSsiIiKS0ZSsiIiISEZTsiIiIiIZ\nTcmKiIiIZDQlKyIiIpLRlKyIiIhIRlOyIiIiIhnt/wPUFO4aRHXRiwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0xb4c3be0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "t=np.linspace(0,0.5,100)\n",
    "C0=[cCr_0,cPy_0,cCrPy_0]\n",
    "C=odeint(Reaction_rate,C0,t,args=(kf_20, pH))\n",
    "\n",
    "plt.plot(t,C[:,0],'r-',linewidth=2)\n",
    "plt.plot(t,C[:,1],'b-',linewidth=2)\n",
    "plt.plot(t,C[:,2],'g-',linewidth=2)\n",
    "\n",
    "plt.title('Cr+Py->Cr-Py\\nk_f=%s  pH=%s'%(kf_20,pH))\n",
    "plt.xlabel('Time (hours)',fontsize=12)\n",
    "plt.ylabel('Concentration (mol/L)',fontsize=12)\n",
    "plt.legend(['Cr','Polymer','Cross-linked CrPy'],loc='best')\n",
    "plt.tick_params(axis='both', which='major', labelsize=12)\n",
    "plt.savefig('Cr-Po.png', dpi=300)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "kernelspec": {
   "display_name": "Python [default]",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.5.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}
